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Abstract 

We analyze correlation functions in a toy model of a random geom- 
etry interacting with matter. We show that in general the connected 
correlator will contain a long-range scaling part. This result sup- 
ports the previously conjectured general form of correlation functions 
on random geometries. We discuss the interplay between matter and 
geometry and the role of the symmetry in the matter sector. 

Introduction 

In theories with fluctuating (quantized) geometry it is non-trivial to define 
the concept of a connected correlation function. This is especially apparent 
in the approaches like dynamical triangulation (DT) which are formulated 
in a non-perturbative, coordinate-free way and permit to go beyond the flat 
background expansion. In such formulations the distance itself is a dynam- 
ical and a highly non-local object. This leads to conceptual and technical 
difficulties when trying to generalize the standard fixed-geometry definitions 

Due to the non-locality the correlation functions are very difficult to 
study and very few analytic results exist |P, Q . A growing body of numerical 
evidence suggests however the existence of a simple structure common to all 
correlators in various ensembles of random geometries . This structure can 
be studied with toy models. 

In this contribution we analyze a simple model of random geometry in- 
teracting with matter which permits the detailed analytical analysis of the 
connected and disconnected correlators. 

The structure of this paper is as follows : we first introduce the model, 
then we derive its thermodynamical properties and general structure of the 
correlation functions in grand-canonical and canonical ensembles. As an ap- 
plication of the developed formalism we then investigate the case of Ising 
spins in magnetic field. We compare our results with the Monte-Carlo sim- 
ulations. 



Branched Polymers 



Branched Polymers (BP) provide a very simple but non-trivial example of 
a random geometry ensemble f^. Moreover they exhibit a wide range of 
properties in common with higher dimensional DT systems 0. 

By Branched Polymers we understand the ensemble of planar, labeled 
trees (T) |^ |^ . Each tree T is weighted with a factor p(T, X) depending on 
the geometry and possibly on the additional "matter" fields X living on the 
vertices of the tree. Partition function of the model is thus defined to be : 

Ter X 

where n = n{T) denotes the number of vertices in the given tree. Considering 
the most general form of two-point actions for geometry and matter fields 
we take p to be : 

p(T,X) = Wvq^,x, n 9m^^^-^q„xy (2) 

i <i,j> 

Here qi and Xi denote respectively the number of branches emerging from the 
vertex i and the value of the matter field(s) in this vertex, <i,j> denotes the 
set of all nearest neighbors pairs. Non-negative coefficients Pg^^ and gq^x;q',x' 
are parameters of the model. The only restriction is that we assume gq,x;q',x' 
to be symmetric ie gq,x;q',x' = gq',x';q,x- The canonical partition function Zn 
is related to ([l|) through the discrete Laplace transform : 

oo 

Z{f^) = Y,e-'^^Zr,. (3) 

n=l 

All the properties of this ensemble can be obtained from the partition 
function of the ensemble Tpi of rooted, planted trees ie the trees with one 
point (root) marked with the condition that this marked vertex has only one 
branch ||^, |] . This vertex thus constitutes a "handle" by which planted trees 
can be glued together. It is convenient to assign to this vertex a dummy 
number of branches go and matter field Xq. In this way we obtain a whole 
set of partition functions for the planted ensemble : 

n 

Zqo,xo{l^)= Yl SZ^"''"!!^''.,^, n ^9.,x,;<7,,x,.. (4) 
TGTp, X i=l <i,j> 

Here n denotes the number of vertices not including the root. The sum over 
all neighbors includes the root (denoted by index zero). 
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Figure 1: Recursive construction of the partition function Zq^^ifJ') of the 
ensemble of rooted planted trees. 



Critical behavior of the model is characterized by singularities of the 
functions Zg^^^^). To find them we first observe that those functions fulfill a 
set of equations (see figure |ip : 



Zq,xif^) = e gg,x;g',x'Pq',x' Z^, J{fi). 



(5) 



Expanding Zg^^^fi) around fi : 
and inserting into (|^) we obtain 

^{5q,x;q'x' " Mg^^.g, y (fl)) Z'g, y (fx) = -Zg^^{fx) 
q',x' 

where operator M is defined by : 

Mg^a:;q',x'{f^) = e"^5'g,x;g',x'Pg',x' " 1 )^«' x'^ (/^) • 



(6) 



(7) 



(8) 



Linear system has an unique solution provided that the operator 1 — M(/i) 
is invertible. The condition 



M(yUo) is i^ot invertible 



(9) 



defines the critical point fiQ and is equivalent to the statement that operator 
M(/xo) has an eigenvalue A*^"-* = 1. As we will see later, this corresponds to 
an elongated or tree phase ||]. 

Around point /iq expansion (j^) is not valid. Going beyond the linear 
approximation, the first non-vanishing term will be in general a quadratic 
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form[| leading to : 

Zq,x{fJ'0 + A/i) ^ ^g,x(/io) - AZg^x{^J'o)^/^■ (10) 



As we show in the appendix ^ AZg ,^.(/io) is the eigenvector associated with 
the eigenvalue A*^°^ = 1. We have chosen the sign in (p!0D as to make the 
AZq^2.(/io) positive. This is ensured by the fact that Zg^^^fi) decreases as /i is 
increased as seen from (1). 

Condition (P) is not the only way the singularity of the partition function 
can arise. The function of Zg^^lfi) defined by possibly infinite series on the 
right-hand side of (^) can develop singularities. This singularities will depend 
on the asymptotic behavior of p's. This scenario corresponds to crumpled or 
bush phase. When both conditions are met simultaneously ie X^^^ = 1 and 
right-hand side of (^ is singular we may have a marginal phase [§, |^ . 

To calculate expectation values of a given operator A we introduce an- 
other kind of a partition function : 

^^(/^) = E E ^) E (11) 

TeT X ieT 

which can be viewed as the partition function of the ensemble of trees with 
one point (root) marked, with operator A inserted at the root : 

G'^(-")= E E^~'>(^'^)^^o,.o- (12) 

We can express G^{ii) by the functions Zg^xifJ') (see figure |^ and reference 
for the relation between planted and rooted ensembles) : 

oo ^ 

^^(/^) = EE^^.--P'^.-^.%(/^)- (13) 

X q=l 

Assuming that we are in the elongated phase we can expand the above 
expression using ([T0|) to obtain : 



1, AZg,^(/i0) 



G^(/^) ^ E E Ag,xP,,xZiM-{i - g ^''77 VA/x) 



9=1 

E E ^5,=^P'/,^^q%(/^o)- exp(-g^^^^^^ v^)- 



(14) 



^Allowing negative values for p's it is possible to make also higher term of the expansion 
vanish and by this tune the system to the multicritical point where the corrections depend 
on A^ii pO|. 
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Figure 2: Rooted ensemble. 
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Figure 3: Correlation function 



Performing the inverse Laplace transform we obtain the large n behavior of 
the canonical function : 

^ _^,on J2 5^ A,,^p,,^ZlJif^o)AZ,Af^o)- (15) 

X q=l 

We get the expectation value of operator (A)„ = ^{J2i^qi,x^)n in canonical 
ensemble by normalizing the function : 

(A) = ^ Ex E^l ^q,xPq,xZ^g~\fXo)AZg.^{fXo) ^^^^ 

Gi ExEJliP5,x^g,x^(/^o)A^9,x(/^o) 

Correlation functions 

We define the integrated (unnormalized) correlation functions of operators 
A and B as : 

TeT X i,j&T 

The distance d{i,j) is the geodesic distance ie the length of the shortest path 
between points i and j. In our case this definition is especially simple as only 
one path exists between two points on a tree. 

To calculate the functions ([17|) we consider the chain of r + 1 points along 
the path joining points i and j and we sum over all possible values of the 
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coordination numbers and matter fields values in each vertex (see figure 
This can be done by the transfer matrix technique. In addition to one-point 
weight each vertex inside the chain with coordination number q has q — 2 
trees attached to it. Those trees can be arranged in g — 1 ways relatively to 



the chain. This gives an additional {q — 1)Z^~ (n) factor per vertex [llO|, [TT 
If in constructing the transition matrix we assign the one-point weights to 
the right-hand side, we end up with a matrix that is identical to the operator 
M defined by except that the first row and column are zero as vertices 
with order one cannot appear inside the chain (see figure We will call 
this truncated matrix M. 



In terms of matrix M the correlation function {\i7\} can be expressed as 
G'^^{r) = yyZ,,{A)(MM^~'M) (18 

a.x a'.x' W - ^)^q\x \N 



with 



Z,,.(A) = e-^'^M,,.p,,. Z«-^(/i). (19) 



The functions Zq ,j,(A) account for the end-point effects and the apparent 
asymmetry is compensated by the asymmetry of matrix the M. The above 
expression can be expanded in the basis of the eigenvectors of the matrix M 
(see appendix ^ for details) : 

(-) = E ^-(A) <i E ^.'-'(B) ^x'- (20) 

k q,x q' ,x' 

The vectors v'^'^^ are the right eigenvectors of the matrix M(/i). They nor- 
malization is defined through the relations (^) and (|60D in appendix 

The formula (^) is exact and valid for any BP system with weights 
defined by (0) . We will now expand this expression around /xq assuming that 
we are in elongated phase. From ( 1T0|) we expect : 

q,x 

As we have shown in the previous section, in elongated phase, Ag°'' = 1 
and oc AZq^xil^o)- Comparing (0) with (|l^) and (|l^) we find that 
A^Q^ = Iq^-* (A) (by 1 we denote a constant operator equal to unity : Iq^^ = !)• 
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Inserting ( |2lD into ( |20D and using approximation (1— x)*" ~ e""^'*" we finally 
obtain the structure of an arbitrary correlation function in the elongated 
phase in the limit A/i —>■ with y/Afi ■ r = const < 1 : 



A;>0 



where 



«^'^ = ^' and 5f = 4y^ (23) 

At this stage we want to make two remarks : i) The first term in (^2|) is a 
function of only one variable x = +6^^)y/Aij, and so the distance scale 

is set by :^7^- As we approach the critical value /xq (thermodynamic limit) the 
size of the system becomes infinite. This justifies the name elongated phase. 
In this phase Hausdorff dimension dn = '2- (this will be more apparent in the 
canonical formulation). In thermodynamic limit shifts can be neglected 
but they are necessary to maintain the scaling at any finite volume. This 



fact was first observed in MC simulations of 2D random surfaces |12]. 

ii) If there exists a finite gap between = 1 and the "next" eigenvalue 
A = lim sup;j>o '^0^'* then the first term will dominate at large distances : 

G'^^ir)^ {A}{B)G^^{r + 6''jl^- + 6^^-) for ^ > ^ = (24) 

with 

6'/^- = 5f ~ 5^^\ (25) 

The function G'-'--'-(r) = G{r) contains no operator insertions. Properly nor- 
malized it gives average volume of the spherical shell of radius r. As such 
it encodes the general information about the average shape and size of the 
system. The formula (|2^ ) states that the whole effect of inserting operators 
is limited to multiplicative factors and shifts of the argument. The shifts are 
additive which means they depend on each operator separately. 

Obviously the formulas (^) and (^) are valid only when A^q\ X^^^ and 

(k) 

X\ are non-zero. We will assume that this is the case for A's and discuss 
the case of A^q^ = 0. Expanding (|2l| ) to the next order : 



z.AM'l ^ -A'^ + A'^^f^ (26) 



q,x 
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we find the corresponding term in the sum (p2D to be : 
W./aTT pvt. -.W 



A-VAm exp ^-a^'^) (r + ^-L) VA^I) = 

We will use this result in the next section when discussing the connected 
correlators. 



Connected correlators 

To study the connected correlators we switch now to the canonical ensem- 
ble. Connected correlation functions G^^(r) are defined in the same way as 
grand-canonical ones, but now summation in (|T^ is restricted to the canon- 
ical ensemble. Additionally we now normalize correlation functions so that 
Ylr^n^i^) °^ ^- Taking the inverse Laplace transform off (|2^) we obtain : 



, r(0) , r(0) 

(28) 



k=i 



where g{x) = xe~^' . The structure of the correlation functions is thus the 
same as for the grand-canonical ensemble with n playing the role of 1/A/i 
and g{x) substituted for e~^. In particular the asymptotic form ( P^ ) is also 
valid in this ensemble. 

We define the connected correlator hyf\ : 

QABconn ^ ^JS^conn-^conn ^ ^( A- ( A) ) (B- (B) ) 

= G^^(r) - (A)G^^(r) - (B)G^^(r) + {A){B)Gl\r). ^^^^ 
The operator = Aq^^ — (A) has a zero average. It is easy to check that : 

.ik)conn _ A^) _ / \ \ -i (k) 



^See and ||] for discussion of other possible definitions. 
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and thence from formula ( pTj) we obtain : 



(31) 

The dots stand for the shifts which depend on A^2^ and B^^ (see (J27\)). The 
shifts '^A{B) obtained by inserting ( pOD into (p3|). 

A more intuitive way of obtaining this resuhs is to insert ( p2l) into the 
expanded expression for the connected correlator (p9|) 0, |13 . 

Comparing (|3T|) with (^81) we see that the scaling long-range part is sup- 
pressed by factor of n. We will call this weak long-range correlations^ The 
non-scaling terms remain of the leading-order in n and for any finite r will 
eventually dominate in the thermodynamic limit. Hoverer when considered 
as functions of the scaling variable x = rj^fn they will vanish in this limit. 

The weak long-range correlations and can be further suppressed by using 
an "improved" version of the connected correlator : 

G^^-"\r) = G^^(r) - {Ay^G'^{r + 5^^ ) - (B)G^^(r + 5^^^') 

+ {A-){B)G'^{r + 5if-+5^4f-). (32) 



Inserting this into (p9|) and omitting the shifts we obtain : 

+ 2v^El«(C-<A>A^)(C-(B>A^)(A<'=))'s(aW^). 

(33) 



We see that the scaling term is further suppressed by a factor ^/n while the 
non-scaling terms remain of the leading-order. 



Ising model 



In this section we consider the Ising model on BP as defined in |T5[. This 
corresponds to the choice of gq,s;q's' = e^"*'* and Pq^s = Pq ■ e~'**. Spins take 

^This is in analogy with the model of spins interacting with infinite range forces but 
with strength diminishing with the system size and leading to similar behavior of the 
connected correlation functions G'™""(r) ~ Q. 
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values ±1. Because we have no geometric interactions, partition function 
Zg^xifJ') does not depend on q. Summing the right-hand side of (|^) over q' 
we obtain [ITH] : 



ZM=e-^(^e^e^F{ZM)+e~'e~^F{Z4^,))^ 



(34) 



where : 



FiZ) = J2piQ)Z''-' (35) 

q=l 

The critical point condition (|^) is now : 

1 = e-^«e^(^e'^F'(Z+(/io)) +e-'^F'(Z_(/io))) 

- e-2'^°2sinh(2/3)F'(Z+(/io))F'(Z_(/io)). (36) 

The transfer matrix M (there is no distinction between M and M as we no 
longer have the dependency on q) is : 

/e^e^F'(Z+) e-^F'(Z_) \ 
M = e"^ (37) 
\e-^F'{Z+) e^e-f^F'iZ^)) 

From (pi]) we have : 
G^^— (r) ^ ^X^B/(a(°)^^) 

+ 2v^FAB(A«)'"^(a«^^) (38) 

We will refer to the first term of the above expression as to "scaling" term 
and to the other as to "non-scaling" one. The expressions for X(y)AB can be 
readout from (pT]). The case of the "curvature-curvature" correlator^ G^ir^ 
is special (operator q is defined by q^^^ = u). We can obtain it directly from 
G^{r) using the relations : 

GAq(O) = G^i(l) and (q) = 2 ^ ^ 



^In DT the curvature is related to the number of ZJ-simphces incident on a (Z? — 2)- 
simplex 
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that are satisfied on any tree. Inserting the above relations into (^) we 
obtain : 



^^"(r) = G^^{r + 2) - 2G^\r + 1) + G^^{r) ^ G"^\r + 1) (40) 

When h = the symmetry imphes Zj^^ji) = Z^{fj,) and the equations (H) 
decouple. The only effect of the presence of Ising spins is to renormalize the 
chemical potential |T^. This leads to : 



Gl^ir) = G~{r) = {tanh (^Y Gl^r) (41) 

where 

g^\r) ^ Cv^re-"^ (42) 

is the point-point correlation function of the pure (without spins) BP model. 
The correlation functions of the operators depending only on the geometry 
remain unaffected by the presence of spins. We see that in this special case 
matter and geometry sectors do not interact : 

Xss = and Fqq = 0. (43) 

The vanishing of Xgs is actually a general feature : whenever symmetry 
enforces the condition Zg^^if^) = Zq{fi) the matter-matter correlators will 
not contain the long-range scaling part : X^x = (see Appendix |C| for the 
proof). 

When h ^ the symmetry is explicitly broken and we have to investigate 
the full system (0,^). Those equations cannot be solved analytically even 
for the simplest choices of p's. They are however not difficult to solve nu- 
merically. Once /io and Z±{no) are known one can calculate other necessary 
quantities. To obtain the i/A/i order terms the standard perturbation theory 
can be used. 

For simplicity we have chosen the weights Pi = P2 = P3 = ^ and all others 
equal to zero. The inverse temperature and magnetic field were set to (3 = 1.0 
and h = 0.1. For those values we obtained Aq ~ 0.5709 corresponding to 
spin-spin correlation length ( 1.785. For a^^-* and and a^^^ we obtained 
respectively 1.1795 and 0.5760. We have compared those theoretical large n 
limit predictions with the results of MC simulations of the systems with 250, 
1000 and 4000 vertices. 

First we have tested the validity of the approximation (0). For r 3> C 
we expect 

G^%r) ^ (s)G(r + Sl^^ ) and G^%r) ^ {s)G%r + 51^^-) (44) 
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n 

range 








theor. 


0.6980 


-0.6533 




250 
10-80 


0.6961(3) 


-0.647(4) 


0.20e-3 


0.6962(3) 


-0.625(4) 


0.90e-4 


1000 
10-160 


0.6979(1) 


-0.644(1) 


O.lOe-2 


0.6979(1) 


-0.652(1) 


0.15e-3 


4000 
10-200 


0.69798(7) 


-0.649(1) 


0.56e-4 


0.69799(7) 


-0.652(1) 


0.26e-3 



Table 1: Fitted values of shifts and average spin values. 

In the table |I| we list the measured shifts and mean magnetization obtained 
by fitting the formula above to the measured values of G^^(r) and Gf^{r). 
The fitting was done by interpolating the functions by cubic splines and then 
minimizing the error function : 

E = - (A)G^^(r - <51^^-))' (45) 



with and (A) as free parameters. As we can see the formulas (|4^) are 
very well satisfied already for the small system sizes. 
For the connected spin-spin correlator we obtain : 

Ass = 0.1414 and Fss = 0.555 . (46) 

The MC results for this correlator are plotted in the figure ^ As expected at 
short distances it is dominated by non-scaling correlations. The line denotes 



the non-scaling part of (|38D in the limit of n — >^ oo. In the inlay we have 
emphasized the scaling part by plotting the correlators scaled by 2y/n as the 
function of the scaling variable x = The line denotes the scaling part 



For connected curvature-curvature correlator we obtain : 

Aqq = 0.6861 and Fqq = 0.0070 . (47) 

We see that in this case the non-scaling part is strongly suppressed. This can 
be seen on the figure ^ where we have again shown the rescaled functions. 
The non-scaling part can be seen only by looking at small r and large n. 
In the inlay we show the small r region. The continuous line denotes the 
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Figure 4: Spin-spin correlation function. In the inlay the rescaled func- 
tions \/nG^{r) are plotted as the function of the scaling variable x = 
The lines denote the theoretical predictions respectively for non-scaling and 
scaling parts. 



non-scaling part of the expression (^) in the n — oo limit. As we increase 
the size we clearly see the emerging signal of the short range correlations 
which approaches the theoretical curve. 



Matter— geometry interactions 

At the first sight the appearance of the long-range component in correlation 
functions of non-critical fields is surprising. It is however clear that it is 
mediated through the geometry. To see this more explicitly we consider an 
example of fields non-interacting with themselves : gq,x;q',x' = 1 but inter- 
acting with the geometry through the term Pg^^- This example can be easily 
solved by introducing new weights and operators : 

P, = J2px,q and Aq = ^^^. (48) 
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Figure 5: "Curvature-curvature" correlation function. In the main part 
rescaled functions ^/nG'^'^{r) are plotted as the function of the scaling vari- 
able X = The inlay shows unmodified functions. The lines denote the 
theoretical predictions respectively for scaling and non-scaling parts 



The correlation functions of the operators Aq with new weight p are identical 
to the correlation functions of the operator with old weights p. We see that 
in this case matter operators behave as geometric operators and although 
matter is non-interacting they pick up long range correlations through cou- 
pling to geometry. 

To see if the effects described in previous section can be explained by 
the above mechanism we assume that the spins are independent and interact 
only with the geometry. The effective coupling term pq s can be read from 

m ■■ 

The effective shift is then given by : 

leff _ 1 (q) (sq)-(q)(s) 
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Figure 6: Comparison of interacting (continuous line) and non-interacting 
(dotted line) spins. 



On the figure ^ we have plotted for comparison]^ (s) ■ and (s) ■ 6s. While 
we see a clear correlation, there is no quantitative agreement. The real effect 
is much more pronounced. This is due to the non-zero spin-spin correla- 
tion length. As a results we have an area-area like interactions rather then 
point-point like as we assumed. When h is increased the correlation length 
decreases and the two values begin to agree as expected (see figure H). Prob- 
ably the proper way to proceed in this case is to use renormalization group 
analysis. Blocking the spins would eventually leave us with a non-interacting 
model and one could expect that the shift would be then given by the formula 
(^). It is so far unclear how to perform this blocking on the BP geometry. 
However a prescription that works for Ising model on 2D random surfaces 



could provide a testing ground for this hypothesis [16 



^(s) • Sg « 5^^^' where (j) — {s + l)/2. This corresponds to the definition used in [Q. 
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Discussion 



We have analyzed in detail a simple model of random geometry and derived 
the general structure of the correlation functions. We observe that in general 
every connected correlator contains a weak-long-range scaling term. This 
term is present even in the correlators of non-critical matter fields. The ap- 



pearance of this term is a direct consequence of the asymptotic behavior (^4]) 
of the disconnected correlation functions. This behavior has a geometric ori- 
gin and is picked up by the matter fields that effectively couple locally to the 
geometry. As we have shown in the appendix P this coupling is suppressed 
in the presence of the symmetry. This is intuitively clear : When h = 
the probability of the spin being up or down must be equal independently 
of the surrounding geometry. When h probability of spin being up is 
dependent on the vertex number q. This can be seen for example from the 
fact that in this case (sq) 7^ (s)(q). 

The long-range scaling term can be strongly suppressed by the use of a 
slightly modified definition of connected correlator (^). This form of im- 
proved correlator postulated in Q is here supported by analytical calcula- 
tions. 

Conversely also purely geometric correlators pick up the non-scaling be- 
havior coming from the matter-matter interactions. In case of the model 
studied here this effect is very small and for non-critical matter it is irrelevant 
in scaling limit. Nevertheless we consider this as an interesting mechanism 
that could play a role in some more complex models which would allow for 
the existence of the critical matter fields. 

While BP provide only a toy model of random geometry, numerical sim- 



ulations indicate that the asymptotic behavior (24) is valid in almost every 
kind of random geometry ensembles. In particular the behavior of the Ising 
model in a magnetic field coupled to the 2D simplicial quantum gravity ex- 
hibits qualitatively the same features and the improved correlator can be 
successfully used Similar structure is also observed for the correlators 
of the local action density of the Abelian gauge fields in the 4D simplicial 
quantum gravity . 



It thus seems that there exists a high degree of "universality" in the 
structure of correlators in systems with a random geometry. Simple models 
as the one described above give us so far unique opportunity to study this 
structure analytically. 
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A Critical point 

Assuming the expansion : 

and inserting it into we obtain following equations : 

J2 (^'/.-;«'.-' - AZ^j^, = 0, (52) 

q' ,x' 

Zq^x = " {^q,x;q',x' " -^(j,x;<j',a;') AZ^, 

q' ,x' 

+ \Y.gq,x;q',x'Pq',x',i(^ " l)(g' - 2) Z^^"? ( AZ^^),, )' • (53) 

q' ,x' 

Equation ( p^ states that AZg)i is an eigenvector of the matrix M(/xo) asso- 
ciated with the eigenvalue A*^°^ = 1. We thus have AZg|2 = Cvf}x- Denoting 
by Vq'^x the left eigenvector of matrix M(/i) associated with the eigenvalue 
A(°) = 1 we obtain from ( [53| ) : 

q,x q,x q' ,x' 

gq,x;q',x'Pq',x'{q' - 1) - Zg^J (v'^^yf ■ (54) 

From this we can easily calculate constant C. 

B Transfer matrix 

Matrix M is not symmetric and it's eigenvectors are not orthogonal. It is 
more convenient to work with a symmetric matrix. To this end we define the 
transfer matrix T by : 

Tq,xW,x'{lJ^) = ^^Mg,3,.gv(/i) q,q' > I (55) 

Jq',x' 
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and zero otherwise, where 



fq,x = Y (g - l)Pg,xZ^/{fi). (56) 

The matrix T(/i) is symmetric, has the same set of eigenvalues A*^*'-' as matrix 
M(/i) and its eigenvectors u'-*^-' are related to the eigenvectors v'^'^^ of the 
matrix M(/i) by : 

vi':l = \^<% ?>1- (57) 

J q,x 



In terms of matrix T correlation functions can be rewritten as : 

Gi'^ir) = E E E E 

q,x T>l,y q' ,x' r'>l,y' 

%^M,.^,, (^'^"')..;.'y iM^)r',y';q',x'^^^. (58) 

Jr,y ' ' ' Jr' ,y' 

This can be further written in terms of an orthonormal basis of eigenvectors 

q,x ■ 



U 



u 



(fc) 



r,y 



q,x;r,y n 
g,x r>l,y J i^^V 

,x'\r',y' J, 

q',x' r'>l,y' •'^''y' 



(59) 



Using (|57|) and the relation between eigenvectors v of matrix M and M : 

vi'l q > 1 

= 1 1 , , ^a. ^ (60) 



Y J2 Mi,x-q',x'vll:l, q = 1 

. ^ q'>l,x' 

we finally obtain the expression 



C Symmetry 

Let us assume that matter field takes values in a group and that the weights 
are invariant under the action of this group : 

9q,y-x;q' ,y-x' 9q,x;q',x' and Pg,y-x Pq,x- (61) 
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From the above and the definition (|I]) we immediately obtain that Zg^x{fi) = 
and so also AZgfi(yUo) = ^^^^i (/^o)- As we have shown in appendix ^ 



this implies Vq^l = v^^^ 



If an operator Ax depends only on the values of the matter fields then : 
J2 ZUA. - {A))vfl = J2 ZqMP.vfl Y.^Ax - (A)) (62) 



q,x 

But from ([T6|) we have 



EgP9,l^g,T^(/^0)AZ,,l(/i0) Ex 1 Ex 1 

which inserted into ( |62D immediately leads to : 

^Z,,.(A-(A))t;W^O. (64) 

q,x 

This ensures that the long-range term in (pOD corresponding to the eigenvalue 
A(°) = 1 will be absent for any correlation function of the type 



r . 
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